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Abstract 

The mathematical theory of quantum feedback networks has recently 
been developed [5] for general open quantum dynamical systems interact- 
^'i . ing with bosonic input fields. In this article we show, for the special case of 

linear dynamical systems Markovian systems with instantaneous feedback 
connections, that the transfer functions can be deduced and agree with 
r^ . the algebraic rules obtained in the nonlinear case. Using these rules, we 

^ l' derive the the transfer functions for linear quantum systems in series, in 

cascade, and in feedback arrangements mediated by beam splitter devices. 
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m . 1 Introduction 

m 

^^ , The aim of this paper is to deduce the algebraic rules for determining the dy- 

f^^- ' namical charactersitics of a prescribed network consisting of specified quantum 

^^ , oscillator systems connected by input-output fields [I], [2]. Physical models in- 

cluded cavity systems or local quantum oscillators with a quantum optical field. 
The resulting dynamics is linear, and the analysis is carried out using transfer 
function techniques [3], [J. The rules have been recently deduced in [5] in the 
general setting for nonlinear quantum dynamical systems by first constructing a 
^ • network Hamiltonian and transfcring to the interaction picture with respect to 

^ I the free flow of the fields around the network channels. However it is of interest 

to restrict to linear systems for two main reasons. Firstly, the derivation here 
for linear systems procedes by an alternative method to the general nonlinear 
case, and we are able to confirm the restriction of the nonlinear formula to linear 
systems yields the same result. Secondly, linear systems are the most widely 
studied models in both classical and quantum dynamical systems theory and 
so it is natural to develop these further. There has been recent interest in the 
development of coherent, or fully quantum control for linear systems [6] -[10] and 
this paper contributes by establishing the algebraic rules for building networks 
of such devices. 
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2 Linear Quantum Markov Models 

The dynamical evolution of a quantum system is determined by a family of 
unitaries {V {t, s) : t > s} satisfying the propagation law V (t3,t2) V (^2,^1) = 
V (ts, ti) where is > ^2 > ^i- The evolution of a state from time s to a later time 
t being then given by ip (t) — V {t, s) 'tp (s). In a Markov model we factor the 
underlying Hilbert space d& \) ® £ representing the system and its environment 
respectively and the unitary V (t, s) couples the system specifically with the 
degrees of freedom of the environment acting between times s and t. For a 
bosonic environment, we introduce input processes bi (t) for i = 1, ■ • • ,n with 
the canonical commutation relations, Jj, 

h{t),bl{s)]=6,,S{t-s). (1) 

It is convenient to assemble these into the following column vectors of length n 

( &i(l) 



b'" [t) 



(2) 



A Markov evolution can be described equivalently by the chronological-ordered 
and Wick-ordered expressions 



rt pt 



V {t, s) — Texp —iT (r) dr = : exp —i / Twick (t) dr : 
where the stochastic Hamiltonian is (with Ej, — Eji and K^ = K) 

n n n 

T (i) = J2 ^v <» 4 (t) hj (i) + ^ F, (^ b\ {t) +Y,^1® ^3 {t)+K<» 1, 

ij = l 1=1 j = l 

and the Wick-ordered generator is given by [13] 

n n 

ijj — 1 i—1 



- Y, Lis,, (g> b, (i) - 2 E ^^^ - '^ 

j=l \ i=l I 



The Wick-ordered coefficients are given by the Stratonovich-Ito conversion for- 
mulae, see appendix, 

1 i. 7? 1 1 1 

S= ?— , L = -i ^F, iJ = £:oo + -ImF ^F^. (3) 

1 + §F' 1 + 5^ 2 I + 5F ^ ' 

Note that H is selfadjoint, and that S* is a unitary matrix whose entries are 
operators on \)'^Y^=\ ^^^Sj^ = % = Y^=\ ^li^kj- In fact, we may write S = 

rp 

e~'"^ with J = 2arctan — . 
2 



In differential form we have 
-y(i,5) = -^ : Twick W V^ (i, s) 



Y, 4 (t) (S,, ~ 6,,)V (t, s) b, (t) + J2 b] (t) UV (t, s) 



i,j=i 



n / n \ 

^ LlS.,V (i, s) b, (0 - 2 E ^1^' - ^^ p it, s) . 

i=i \ i=l / 



Note that all the creators appear on the left and all annihilators on the right. 
This equation can be interpreted as a quantum stochastic differential equation 

[I], [m, [El- 

We sketch the system plus field as a two port device having an input and an 
output port. 

system 
output, b°"' I 1 input, b'" 



Figure 1: input-output component 

The output fields are defined by 6™* {t) = V (t, 0)'^ b, (t) V {t, 0) and we have 
the input-output relation 

n 

where Sij (t) = V (t, 0)^ S,jV {t, 0) and Li (t) ^ V {t, 0)^ L,V (t, 0). More com- 
pactly, b™' (t) = 5 (t) b'" {t) + L{t). 

Let X be a fixed operator of the system and set X {t^to) ~V (t, to) XV (t, to), 
then we obtain the Heisenberg-Langevin equation 

^X{tM) - V{tM)''-[X,T{t)]V{tM) 
at I 



= b] (t) V (t, to)^ [SlxSk, - 5^,X) V (t, to) b, (t) 

+bl{t)V{t,to)^Sl[X,Lk]V{t,to) 
+V {t,to)HLlx]SvV it,to)b, (t) 

+V{t,toyihl[X,Lk] + ^[Llx]Lk-i[X,H]\v{t,to). 

Note that the final term does not involve the input noises, and that the expres- 
sion in braces is a Lindbladian. In the special case where S = 1, this equation 
reduces to the class of Heisenberg-Langevin equations introduced by Gardiner 



2.1 Linear Models 

We consider a quantum mechanical system consisting of a family of harmonic 
oscillators {oj : j = 1, • • • ,m} with canonical commutation relations [0^,0^] = 







al'4 



and 



aj,al 



— djk- We collect into column vectors: 
ai 



(4) 



Our interest is in the general linear open dynamical system and here we 
make several simplifying assumptions: 

1) The Sjk are scalars. 

2) The L' s are linear, i.e., there exist constants Cjk such that Lj = ^j, CjkCik- 

3) H is quadratic, i.e., there exist constants tUjk such that H = X^iTc 'K^jkO.k- 

The complex damping is ^L^L + iH = —a^Aa where A = —^C^C — ifl 
with C = (cjk) and Q — (cojk)- The Hcisenberg-Langcvin equations for a(t) = 
V (t, 0) aV (i, 0) and input-output relations then simplify down to 

a(i) = Aa{t)-C^Sh{t), (5) 

b°"*(i) = Sh{t) + Ca{t). (6) 

These linear equations are amenable to Laplace transform techniques [3] , [4] . We 
define for Re s > 



oc 

St 



C{s)= / e-'''C{t)dt, (7) 







where C is now any of our stochastic processes. Note that a (s) = sa (s) — a. 
We find that 

a(s) - -{sI„,-A)-'C^S9''{s) + {sI,n-A)-'a, 
b°"* (s) = S-b'" (s) + Ca (s) . 

The operator a (s) can be eliminated entirely to give 

b°"'(s)=S(s)b'"(s)+C(s)a (8) 

where the transfer matrix function is 

E{s) = S-CisI„,^Ay^C^S (9) 

andf (s) =C(s/„,. -A)'\ 

As an example, consider a single mode cavity coupling to the input field via 
L = ^/7a, and with Hamiltonian H = coa'^a. This implies K = ^ + iu> and 



C = ^7- If the output picks up an additional phase S 



— pi4> 



transfer function is then computed to be 

^cavity \^ ) 



s + iw + ^ 



the corresponding 



(10) 



2.2 The Transfer Matrix Function 

The models we consider are therefore determined completely by the matrices 
{S,C,n) with S e C"^",C e C"^™ and Vl e C™''™. We shaU use the con- 
vention -^ — ^ (s) = D + Cis-Ay'^B for matrices A £ C™><",B g 



C" 



,Ce 



and D e 



^, and write the transfer matrix function as 



Eis) 



■ A 


~c^s ■ 


C 


s 



is), 



(11) 



where A — ~^C^C — iVt. We note the decomposition 



C {sl„, ^ A)-^ C^ 



In the simplest case of a single cavity mode we have 



A 


-ct 1 


C 


/„ J 



s. 



^-•cavity 



(s) 



V7 



-V7e^^ 



(s). 



Lemma 1 For each w G R, i/ie transfer Junction S (iw) = .^ (0''" + iw) is uni- 
tary whenever it exists. 

Proof. The decomposition follows immediately from (fTTj) . Wc have then 
for instance 



S(0+ + 
= I-C- 



s(o+ 



1 



I-C 



1 



iCtC + ifi' 



ct 



I-C 



1 



iCtC-if7' 



c-t 



, , ictC + if7' + ictC-2f^'-CtC, , , 

iCtC + if7' 12 2 1 iCtC-if}' 



1 



ct, 



where Vl' = fJ + cj. The term in braces however vanishes identically, leaving 
S (0+ + iu}) S (0+ + iuf = /. The relation S (0+ + iw) S (0+ + iu)^ = / is sim- 
ilarly established. ■ 

Whenever appropriate, we may determine S from its (unitary) values on the 
imaginary axis by using the Hilbert transform 



^(^)-i^^ 



S(itj) 
Lo -\- is 



dio. 



In general, the real and imaginary parts of A need not commute - that is, 
[C^C, O] need to be identically zero. However, when this does occur we recover 
a multi-mode version of the cavity situation. 



Lemma 2 If A is a function ofC^C then 

s- A^ 

where A is a function ofCC^ and S may be analytically continued into the whole 
complex plane. 

Proof. Here we must have A = —^C^C— ie {C'^C^ where e is a real valued 
function. We set A = -\CC^ - ie {CC^). From the identity Cf {C'^C) C^ = 
CC^ f (^CC^) for suitable analytic functions /, we have 



s-A]Eis) = {s-A 



I ^CC^ 

s-A 



S= s-A-CC^]S 



however, —A — CC^ = A^, and this gives the result. 

The hermitean matrices C^C and CC^ will have the same set of eigenvalues: 
to sec this, suppose that (/> is a non-zero unit eigenvector of CC'^ with eigenvalue 
7, then ?/; — ^~^/'^C^ ip is a unit eigenvector of C^C with the same eigenvalue, 
conversely, every eigenvector ^/^ of C^C with non-zero eigenvalue 7 gives rise to 
a nonzero eigenvector (p ~ j^^^^Ctp of CC^ . 

Let CC^ have the spectral form ^^ Ik^k with real eigenvalues 7^. and cor- 
responding eigenprojectors Ek, then we have 



1. 

S + i7fe -t- lEk 



Eis)^Y. :^i!' ^^ ^^^^ 



where Sk = s (7^). In particular, the rational fraction is of modulus unity for 
imaginary s (= itu) and we may write 

k 

where 4>f. {uj) = Q-^'S i(^+e'°)+^'°/2 - Note that S (0+ + itu) is clearly unitary and 
the limit w ^ is well-defined. This limit will equal —S in the special case 
that K is selfadjoint (i.e., e = 0). S may be analytically continued into the 
negative-real part of the complex plane. The poles of S then form the resolvent 
set of K, and the zeroes being the complex conjugates. ■ 



3 Introducing Connections 

The situation depicted in the figure below is one where (some of) the output 
channels are fed back into the system as an input. Prior to the connection 
between output port(s) s\ and input port(s) r; being made, we may model the 



component as having the total input b" 






and total output b° 



uout I where the b'" and b°"* may be multi-dimensional noises (we in fact 
only require the multiplicities to agree for j = i,e respectively). 




figure 2: A quantum system with feedback 
The transfer matrix function takes the general form 



A 






Y^jC.Sj\ lli^j^: 



je 






Si, 



When we make the connection, we impose the various constraints ^^c^.-, (t) = 
KT) (* ^ ^) "^here output field labelled s; (j) is to be connected to the input 
field r; (k) where r > is the time delay. We assume the idealized situation of 
instantaneous feedback t ^ 0+ . To avoid having to match up the labels of the 
internal channels, it is more convenient to introduce a fixed labelling and write 

b°"' {t-) ^ Tybi" {t) 

where rj is the adjacency matrix: 



Vsr = 



1, if (s,r) is an internal channel, 
0, otherwise 



The model with the connections is then a reduction of the original and the 
remaining external fields are the input b^" and output b°"*. 

Theorem 3 Let (rj — ^ij) be invertible. The feedback system described above 
has input-output relation b"""* = Srcdb^ + ^^cd^ '^'^'^ ^^^ reduced transfer matrix 
function 



' — ^rcd 



^rcd 


Cj.g^S'rod 


Crcd 


Sicd 



where 



Srcd — '-'ri 



'-'rod — »5ee i »^ei [V ^u) "^ie: 
Cred = Seiiv-S-.y^Q+Ce, 

A,ed = A-Y,Cpf,{v-Sv,)-'Q. 



1 



A, 



red 



(12) 



Proof. The dynamical equations can be written as 

k{t) = Aa{t)-J2c}S,kK{t), 

fe=Le 

Now the constraint rjh™ = b°"* imphes that 

bj" (t) = (77 - S,y' (5ieK" (t) + Qa (t)), 



and so 






J=i,e 



or 



a (s) = ^-— J2 C] is,, + 5,i (ry - ^i,)"! 5ie) be {s) + ^a, 



J=i,e 



with Arod as above. Consequently, 



— Ogj Dj -\- Dee t)e '^ Oea 
= Orcdbe + CfodS 
= ^rcdbg + Crcd^ 



where 



Jrod / ^ L/fcd T C^ ( '-'je + '-'ji V*? ^ "^iij "^ie j ) 



^red 

ircd 

Srcd ^ L/i-cd" 



and S'red, C-ccd are as in the statement of the theorem. 

We now show that Y.j=\^e C][Sje + S,; (ry - Sy,y^ 5ie] = Cl^S.^^- Now 

^CJ[S'je + 5,i(77-5ii)-'5ie] - q^[5ie + 5ii(77-^ii)-^^ie]+Ct^rcd 

= q^r;(ry-^ii)-'5ie + Ct^rcd, 
while C^^d-^rcd - C}{t^^ - SlX^SlS,,^ + Ct^rcd. However, 

(r;t - 5it)-i5it^,ed - (ryt - ^i|)-i5't(5ee + 5ei (r, - ^i,)-^ ^le) 



and using the identities SlS\\ + Sl-Se\ ~ 1, SlS\e + Sl-See = 0, this reduces to 
(Tyt - Sly^SlS,,^ = (r,t - 5.t)-i ^-SlS;, + (1 - SlS,) (77 - S.y' S;, 

= iV^ - sir' [si (fj - 5ii) + (1 - SlSy,)] (rj - Sy,)-' 5ie 

= ri{ri-S]\y S]e- 

Therefore S,.od = ■S'rod - !],=: e^iod 1 — C^S^cd, as required. 

■' ' s- A^cd 

For consistency, we should check that we have A^cd = —■^Cj.^^Crcd — if^rcd 
with rij-od selfadjoint. Indeed, setting A = —^CIQ — ^C'lC'e — i^ and substi- 
tuting in for Crod and K^-^d we find after some algebra that 

f^red = fi + Im {c/5ii (t? - Sii)"^ Ci} + Im {c^S-ei {tj - Sv,y' C,} . 

The manipulation for this is trivial except for the calculation of the term of the 
form iC| XCi where 

X = l + 2Sv,iv-Sy,)-'-{TJ^-Sl)-'SlS,;{v^Sy,r' 



i^-vsi) 



t^-l 



Sy^^ - vSl 



2ilmr — M— r =2ilm 



1 - Sy,rj^) 
{Syiv-Syy'} 






L 






1 - Sy,TJ^ 

where again we use the identity SlSy, + Sl^Se, — 1 
In terms of the parameters {S, L, H) with S 

I „' I and H = a^ila, we have that the feedback system is described by the 
reduced parameters (S'lod, -^red, -f^rcd) where 

"Jrcd = '-'ee + '-'ei (V — by) Oje 
Lied — i>e\ (rj — Syj L\ + Lg, 

H,ed = i/ + Im{Lt5ii(7/-5ii)-iLi}+Im{LtS'ei(?7-5ii)-^Li}. 

(13) 

The same equations have been deduced in the nonlinear case by different argu- 
ments 0. Note the identity Im ^LJSy [t] - Syy' l\ = Im (l^ (?/ - Syy' l\. 



Remark 4 Let U be a unitary operator on a fixed Hilbert space ^ — ^i 

C/11 C/12 

U21 U22 
transform ffy' is the superoperator defined by 



SJ2 which decomposes as U 
s the supei 

iffy' (x) = Uii + c/12 (1 - xu22y' XU21 



The non-commutative Mobius 



defined on the domain of operators X on 9)2 for which the inverse (1 — XU22) 
exists. The transform ifiif^^ maps unitaries on S)2 iri its domain to unitaries in 

Remark 5 In particular, S^cd is unitary as it equals if'cj*^ (^) where £, = rj^^ 
with r/ being unitary. We may expand the geometric series to write 

00 

'-'rod = "Jee + '-'ei^'Jie + i'eiC'-'iiC'-'ie + '-'ei^'-'iiC'-'iif "-"ie + ' ' ' = Oee + / ^ "-"eiC ('-'iiC) '-'ie 

which shows that S^cd can be built up from contributions from the various paths 
through the network. Likewise 

00 

CXj 00 

ff + ^Im{Lt(5i;erii} + 5]lm{it^eie(^iie)"ii}. 



H, 



red 



ri=0 



n=0 



4 Systems in Series 

As a very special case of feedback connections we consider the situation of 
systems in series. This is referred to as feedforward in engineering. 



0^2 ^^O 



■0^1 '^lo 



figure 3: Cascaded systems 
The individual transfer functions before the connection e = (si, r2) is made 



are given by ^^ — 



concatenated to give 



[ A, 


-cjs, ] 


with Ai = 


^c^'c 


[ a 

ve 


S^ \ 


2^i ^« 




' A1+A2 


-cls. 


-ClS2 ' 


S == 


Ci 


Si 







C2 







^2 



CJCi — iQi.and these may be 



To use the formula for the reduced transfer function following connection, 
we must first of all identify the internal (eliminated) and external fields: here 



bi' 



b 
bi' 



? I b° 

111 I ' *-* 



"UOUt 

b°"* 



hf" 



V.out / ' 
"2 



10 



L\ — Li, Le — L2, 



and 

Se, See )~[S2 

with trivially r/ = 1. The reduced transfer function is then readily computed to 
be 



" Ai + A2 - C^2S2Ci 


- {clS2 + cl) Si ' 


C2 + S2C1 


S2SI \ 



Likewise we deduce the relations 

S — S2S1, L ~ L2 + S2L1, 



H = Hi+H2 + Iin{LlS2LiY (14) 



The same equations have been deduced in the nonlinear case by different argu- 
ments [5]. 



4.1 Feedforward: Cascades 

If the two systems are truly distinct systems, that is, if they are different sets 
of oscillators, then we are in the situation of properly cascaded systems. In 
this case one would expect that the transfer function to factor as the ordinary 
matrix product Sgcrics = 22S1. We now show that this is indeed the case. 

Lemma 6 Let Sj be transfer functions for rrij oscillators coupled to n fields 
(j — 1, 2). If we consider the ampliated transfer functions for mi+m2 oscillators 
coupled to n fields 

—C^Si 
^1= \uu/v 




then 



^series 



(15) 



Proof. Wc compute this directly. 

' Ai 

-C2S2C1 —A2 



(Cl,C2) 

= 52^1+ (Ci,C2) 

= S2S1 — (Cl,C2 

= S2-C2 [s - A2 
giving the result. ■ 




C2S2S1 



11 



5 Beam Splitters 

A simple beam splitter is a device performing physical superposition of two 
input fields. It is described by a fixed unitary operator T — I \ G U (2): 



^,out 

V-out 
"2 



M 



bJ," 



This is a canonical transformation and the output fields satisfy the same canon- 
ical commutation relations as the inputs. The action of the beam splitter is 
depicted in the figure below. On the left we have a traditional view of the two 
inputs being split into two output fields. On the right we have our view of 
the beam splitter as being a component with two input ports and two output 
ports: we have sketched some internal detail to emphasize how the scattering 
(superimposing) of inputs how ever we shall usually just draw this as a "black 
box" component in the following. 





Figure 4: Beam-splitter component. 

To emphasize that the beam splitter is an input-output device of exactly the 
for we have been considering up to now, let us state that its transfer matrix 
function is 



^bcam splitter 



■ 








T 



= T. 



Our aim is to describe the effective Markov model for the feedback device 
sketched below where the feedback is implemented by means of a beam splitter. 
Here we have a component system, called the plant, in-loop and we assume that 

it is described by the transfer function Sq 



b™ 



[ Ao 


-cls, 1 


[ Co 


So \ 



hf 



"UOUt 

"2 



/ 


plant 










in 

2 
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Figure 5: Feedback using a beam-splitter. 
It is more convenient to view this as the network sketched below. 

plant 



■0^3 S3o 



b^* 



O 32 r2o. 



^sin^ 



beam splitter 



Figure 6: Network representation. 

Here we have the pair of internal edges (5277*3) and (53, r2). The transfer 
function for the network is 



■-unconn. 



^0 








— Cq5o 





Til 


ri2 








T21 


i 22 





Co 








^0 



with respect to the labels (0,si, 52,53) for the rows and (0,ri,r2,r3) for the 
columns. This time the external fields are bf = b'j", b°"* = b°"* = Tiib'j^ + 
Ti2b2" while the (matched) internal fields are 



That is 



b°"* = 



T22 
So 

^ei-(Ti2,0), 

La 




Kout I — 



_ ( T2lW^+T22ht 

5*0 bjj" + Lo 



S; 



u == 






T21 


7^11, 



ii = 0, 7/ = 



1 

1 



Substituting into our reduction formula we obtain 



T11+ (Ti2,0) 



1 

—5*0 

Til + 2^12 (1 — SctT22) S0T21, 



—T22 

1 



T21 





C = 



-T22 
1 



1 

-So 



{Ti2,0) 
= T12 (1 — 50T22) Co 
O = r^o + Im (0,4) 

= r!o + ImCt (1-^0^22)"' Co 





Co 



— 722 
1 



1 

—5*0 





Lo 
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and so, when the connections are made, the transfer fmatrix function is 



^0 — Cr,SoT22Cn 



Ti2 (1 — 50122) Co 



-CqSoT2i — CqSq (I — S0T22) T22 



Tn + T12 (1 — S'oT22) SqT2i 



Note that S = Lpp*^ (So) where ^|ri (z) = Tn +T12/3 (z^^ - r22) T21 is the 
Mobius transformation in the complex plane associated with T. 

I , and X + iy — Soiy, then 



If we further set T 



C^C 



1 - S^v 






rCnCo 



1 — x^ — y^ 



ImCt (1 - Squ)-^ Co = Im (- ^—} C^Co = - 



Co Co, 

:ClCo. 



(l-x) +?;: 

In particular, if we take a single oscillator in- loop with Sq — e"^" , then we 
obtain S = e*"^ and the phase is determined by the Mobius transformation. 
If we further have Lq — ^/Toa, ^o — ojoa^a, we find that L = e^\/^a and 
H — Lua^a where 



7 



1 - X — y 
(l-xf +y' 



:7o- 



y 



(l-x) +y' 



a P 



-^0, 



with S'o = 1, Wo = 



and (5 is a real phase. In the specific case T — , „ 

\ p —a 

considered by Yanagisawa and Kimura [3], we have x ^ —a and y = 0, therefore 

we find 

1-a 



!'■ 



which agrees with their findings. 



1 + « 



7o, w = 



An alternative computation of S is given by the following argument. We 
consider the input-output relations 



and eliminating b^"' = (1 - 722S)" 



T2ihT + T22eoao 



yields 



br* = \Tii + T12S0 (1 - r22So) ' T21J bi" + T21 (1 - S0T22)"' eoao- 

That is 

S = Tn + T12 (So 1 - T22)"' T21 = </p2ri (So) . 

We remark that if T12 and T21 are invertible, then we may invert the Mobius 
transformation to get 



"-1 

"0 



T22 + 721 



1 



S-Tn 
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To illustrate with a cavity mode in-loop, we take the beam splitter matrix 



to be T 



a /3 



with or 



P' 



1, and the transfer function ^o (s) 



1 + aSo 



1 — g 7 



l-g 7 ■ 
1+a 2 



6 The RedhefFer Star Product 

An important feedback arrangement is shown in the figure below. 

A 

bin ^ c= "uout 



V-OUt 

"3 



V-out 
D4 



-0^2 S20- 



-OS3 rac^ 



-OS4 r4Q, 



b°"* = hf 



B 



Figure 7: Composite System 



described 



, Q,A and B by 



We shall now derive the matrices for this system taking component A to be 

qA qA I ^ I r^A 
'-'21 '-'22 / V "-^2 

J7b- The operators of systems A are asumed to commute with those of B. 
We have two internal channels to eliminate which we can do in sequence, or 
simulataneously. We shall do the latter, here we have 



cB 
■^33 
oB 
"^43 



cB 

"^34 
oB 

O44 






Si. 



qA 



qA 







5*44 



qB 



,Sii 



qA 



qA 

"-"22 





qB 

"-"43 



qB 



and 



U = 






Li = 



-^3 



?? = 



1 

1 
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The parameters are therefore 



S. 



a 



oA 
■^11 





O^ = il^ + fii 



+C^' (1 





qB 



oA 
•^12 





oB 



_qA 
"-"22 



•^33 



oA 
"-"21 





-■34 



c/l I qA qB (-1 qA qB\ 1 cA 
Oil + '-'l2'-'33 V^ '-'22'-'33; '-'21 



^4^3 (1 



cA cB\ 
'-'22'-'33j 



'-'21 



•^12 (1 '^22'5^33) 
544 + 5*43 (1 — 5*22 S'33 



S: 



34 
'599 S'; 



22'-'34 



"-"12 
S. 



43 



Ci^ 

^4^ 



^1 + '5l2'533 (1 '522'5: 

C? 



qA 
'-'22 



^9^ 



33 






'-'12 V '-'22'-': 



33 j 



qB (-1 qA C-B^^l r^A I cB cA ('1 _ qA cB\ 
'-'43 V^ '-'22'-'33; ^2 f '-'43'-'22 l,^ '^22'-'33^ 



Ci 



Im 



oA c'B\ 
'-'22"-'33; 



{^3^^ (1 

^ qA r-iB 
'-'22^3 

r'ci 



qB cAy 



^3^ 



^2 V-^ ^ ^22"^33J 
r<Af qA (-1 qB qA 
^1 '-'12 {^ '-'33'-'22 



St 



(1 



qB qAy 
'-'33*^22J 



qA r^A 
'-'33 '-^2 



c^'{\-s^,siy 



c 



qb r<B 

'-'33<-'2 



r'Al cA (-[ cB q/ 

-L^i O12 {^ '-'33'-'2 

r'^tcB/'i qA qb\^^ qA r-iB , ri^'< qb (-1 qA cb\^^ riA\ 

O4 043(^1 '522'-'33J '^22'-'3 +'-'4 '-'43 1^ '^22'-'33J '-'2 J • 



7 Appendix (Stratonovich to Ito Conversion) 

It is convenient to introduce integrated fields 

B,it)= f h{s)ds,B\{t)= f blis)ds,A,,it)= f blis)b,is)ds. 



Jo Jo Jo 

Bi (t) and Bj (t) are caUed the annihilation and creation process, respectively, 
for the ith field and coUectivey are referred to as a quantum Wiener process. 
Aij (t) is called the gauge process or scattering process from the jth field to the 
ith field. A noncommutative version of the Ito theory of stochastic integration 
with respect to these processes can be built up. The quantum Ito table giving 
the product of infinitesimal increments of these process is 



X 


dBk 


dAki 


dBJ 


dt 


dBi 





5,kdBi 


Sadt 





dh-ij 





SjkdAii 


SjidBJ 





dB] 














dt 















The Ito equation for the unitary process is then dV — (dG) V where 

n n n / ^ n \ 

dc = J2 (^^^ - '5^^)'^^^^ + E ^^^^l - E ^l^^i^= " 2 ^ ^^^' " '^ '^*' 

i.j — l i— 1 i— 1 \ ^—1 / 

The Stratonovich form is dV = —i (dE) o V where 

n 71 n 

dE= J2 E^dA.j + J2 PidBJ + J2 ^j'^^J (*) + ^^*- 

i,j=l 1=1 j=l 
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and we define the Stratonovich differential to be {dX)oY = (dX) Y+^ (dX) (dY) 
with the last term computed using the Ito table. We have the consistency con- 
dition dV = {dG) V=-i {dE) y - § {dE) {dG) V or 



dG ^ -idE - 


- \ {dE) {dG) , 




and using the table we see that 






S-l = 


-iE-'-E{S- 


- 1 


L = 


-^f-\el 




-]-L'^L-iH = 


-iK~'-F^L 





which can be solved to give the relations ([3]). 
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